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(57) Abstract 

A method is disclosed of determining 
position uncertainty of directional boreholes 
A, B in which the true position of well A, r 1 , 
will lie in the ellipsoid centred about rl, and 
the true position of well B, r 2 , will lie within 
the ellipsoid centred around r2. The method 
comprises measuring either the azimuth and 
inclination of the borehole at intervals along 
it, or the total displacement from the surface 
location of successive points, determining the 
magnitude of systematic and random errors 
inherent in the measurements, determining 
the variance matrix describing the position 
uncertainty of each station in the survey, 
determining the point of closest approach of 
the borehole to a point of interest, determining 
the variance matrices representing the position 
uncertainty of each of the point of interest 
and the point of closest approach, determining 
the vector between the point of interest and 
the point of closest approach, determining 
the normalised distribution of the length of 
the vector according to the variance matrices, 
and integrating the normalised distribution to 
obtain the probability that the borehole lies 
within the specified radius of the point of 
interest. 
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A PROCESS FOR DIRECTIONAL DRILLING 

The present invention relates to methods of determining 
the positional uncertainty of a directional borehole, and 
relates particularly , but not exclusively, to methods of 
determining positional uncertainty of a borehole of an oil or 
gas well . 

Controlled directional drilling is ' the process of 
maintaining the trajectory of a borehole along a predefined 
course from a surface location to one or more subsurface target 
locations. Directionally drilled boreholes are widely used for 
the exploitation of subsurface natural resources such as oil 
and gas, and coal bed methane, and are also used during the 
construction of so-called "freeze rings" whilst sinking mine 
shafts through water bearing formations, for scientific 
purposes including the analysis of ground water flow and the 
performance and evaluation of underground nuclear tests, as 
well as for the controlled extraction of ground water 
contaminated with environmentally damaging chemicals and the 
creation of conduits for pipes or cables beneath rivers, lakes 
or other obstacles. 

A variety of techniques may be employed to deviate the 
path of a directionally drilled borehole. When rotary 
drilling, buckling of drill collars within the bottom section 
of drill string (the "bottom hole assembly (BHA) 11 ) can be 
controlled through careful selection of the position and gauge 
of stabilisers, and thus the tendency of the BHA to increase, 
maintain or decrease its angle from the vertical may be 
changed. The stabilisers employed may be of fixed gauge or 
have gauge adjustable from the surface, and in either case 
control over inclination is more readily achieved than 
direction. Downhole motors are commonly used to achieve better 
control over direction, and power the drill bit so that the 
drill string need not be rotated. By introducing a "bent sub" 
into the BHA above the motor it may be orientated with respect 
to the axis of the borehole and thus induced to drill a smooth 
arc . 
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Controlled directional drilling requires that the path of 
the borehole is surveyed during drilling in order that 
appropriate steps may be taken to maintain the well path along 
a desired trajectory. It is extremely difficult precisely to 
control the buckling or orientation of a BHA and thus 
accurately drill in a predetermined direction. Similarly, 
measurements of the borehole vector, i.e. couplets of 
measurements of borehole inclination from vertical and azimuth 
(the observed borehole direction measured from north) are 
subject to errors in inclination, azimuth, and along-hole 
depth. Computation of position along the. borehole thus leads 
to an accumulation of error so that uncertainty in the position 
of stations (i.e. points at which observations of borehole 
vector are made) is greater in stations further along the 
borehole, and lesser in stations closer to the origin. 

Borehole position uncertainty causes considerable 
problems for applications of controlled directional drilling. 
In particular, the economic development of oil and gas fields 
off shore, and in the Arctic and other environmentally 
sensitive areas, frequently calls for a controlled directional 
drilling of multiple wells radiating from single platform, sub- 
sea template or rig site. In such circumstances, wells are 
closely spaced at the surface, and there is a serious risk that 
any new well may intersect an existing well, any such collision 
risking the release of hydrocarbons into the environment, and 
subsequent fire and pollution. The actions required to 
minimise the impact of such a collision include "shutting in" 
the well, i.e. ceasing production or injection and eliminating 
excess pressure within the borehole, and adopting more cautious 
and time consuming drilling methods. Both of these options are 
extremely costly. 

A further problem of borehole position uncertainty is 
encountered when attempting to deliberately intersect an 
existing borehole. Such action may be necessary if control of 
an existing producing well is lost in a "blow out", and has 
also proved necessary when re-entering an abandoned well in 
order to undertake fresh operations. In both cases, lack of 
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knowledge of the borehole position requires that expensive 
"well location" logging equipment be run. 

Computer models have been used for a number of years to 
predict borehole position uncertainty and the results used to 
reduce the risk of well collision by computing the separation 
of boreholes from one another and accounting for their 
positional uncertainties in an empirical manner. Such methods 
have been introduced to systems of control employed at rig 
sites to provide regular "go / no-go" decisions during the 
critical section down to and immediately below the "kick-off". 

Such borehole uncertainty models tend to be variations on 
models first published by Walstrom J E # Brown A A and Harvey R 
P - An analysis of Uncertainty in Directional Surveying - 
Journal of Petroleum Technology, April 1969, PP515-522 or, more 
commonly, Wolff C J M and de Wardt J P - Borehole Position 
Uncertainty, Analysis of Measuring Methods and Derivation of 
Systematic Error Model - Journal of Petroleum Technology, 
December 1981, PP2339-2350. Such methods may be used to 
develop ellipsoids bounding the volume within which the 
borehole location may be expected to lie. 

The primary method of borehole intersection avoidance 
currently used is to generate such ellipsoids or ellipses at 
regular intervals along drilled and proposed boreholes and to 
compare the separation of the boreholes to the size and 
position of the ellipses. This comparison may be undertaken 
visually, or through identification of the point of closest 
approach of the proposed borehole to the drilled borehole and 
subsequent computation of a simple ratio known as the 
separation or clearance factor defined by F = R / (r x + r 2 ) , 
where r 1 and r 2 are the radii of the ellipses of uncertainty on 
the proposed and drilled boreholes respectively along the 
direction between them and R is the separation along that same 
direction. 

Borehole collisions are most likely to occur where 
separation is low, which is generally near the surface. Since 
there the boreholes are nearly vertical and their ellipses of 
uncertainty approximately circular, a factor F greater than 1 
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implies that the two ellipses of uncertainty do not overlap one 
another and if, as in the method of Wolff and de Wardt above, 
the ellipses are held to bound all possible positions of the 
borehole location, then collision of the two boreholes cannot 
occur. At higher borehole inclinations, the ellipses are no 
longer approximately circular, and a separation factor greater 
than 1 no longer implies that the ellipses cannot intersect. 
In such circumstances, r ± and r 2 may be replaced by the major 
axes of the two ellipses to produce a "fail safe" solution. 

Typically, criteria are set down defining actions to be 
taken during drilling at values of F approaching 1 . For 
example, in order to continue drilling at F = 1.5, it may be 
necessary to shut in and bleed off excess pressure in an 
adjacent well. At F = 1.25, surveys of the borehole being 
drilled may be required to be taken more frequently or with 
higher precision sensors. When F is less than 1, it may be 
deemed necessary to cease drilling, plug the borehole and 
deviate it towards a safer area . 

All of the above practises suffer from the disadvantage 
of considerable expense in terms of time, equipment, and 
deferred production from the well. In particular, prior art 
techniques suffer from the disadvantage that they do not permit 
a quantative estimation of the risk of borehole intersection, 
and can therefore cause safety procedures to be prematurely 
invoked such that substantially more cost is incurred than the 
actual risk justifies. 

This inability quantitatively to estimate the risk of 
borehole intersection has substantial cost implications. The 
slot separation on off shore structures and sub-sea templates 
is typically of the order of 10 feet, largely to ensure that 
the ellipses of uncertainty of the boreholes drilled therefrom 
do not intersect. For many off shore structures, especially 
those to be drilled from jack-up drilling units, the topside 
weight is proportional to the area, which is in turn related to 
the square of the slot spacing. Likewise, the jacket weight is 
proportional to the top side weight, and the ability to reduce 
slot separation without compromising safety would be highly 



WO 96/35859 PCT/GB96/01 123 

5 

desirable in terms of the substantial capital savings gained. 

Preferred embodiments of the present invention seek to 
overcome the above disadvantages of the prior art . 

According to an aspect of the present invention, there is 
provided a method of determining position uncertainty of a 
directional borehole, the method comprising: 

measuring azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement from a surface location of 
successive points using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made and their relationship with 
azimuth and inclination of the borehole and displacement from 
the origin and elapsed time at each successive point of 
observation; 

determining a first variance matrix representing position 
uncertainty at each location at which said measurements are 
made ; 

determining the point of closest approach of the borehole 
to a reference point; 

determining second variance matrices representing the 
position uncertainty at said reference point and said point of 
closest approach; 

determining a vector between the reference point and the 
point of closest approach; 

determining a normalised distribution of the length of 
said vector according to the second variance matrices; and 

integrating the normalised distribution between limits of 
zero and a value^ equal to the sum of a predetermined radius and 
the radius of the borehole, wherein the integral obtained 
represents the probability that the borehole lies within the 
predetermined radius of the reference point . 

This provides the advantage of permitting a quantative 
assessment of the risk of borehole intersection to be made 
which in turn enables borehole separation to be reduced without 
compromising safety. 

In a preferred embodiment, the step of determining the 
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first and / or second variance matrix and / or the normalised 
distribution are performed by application of a stochastic 
technique, and preferably by application of Monte Carlo 
modelling, and the integrating step performed by examination of 
the frequencies thus derived. 

In a preferred embodiment, the first variance matrix is 
determined from estimates of expected drilling precision. 

Preferably, the reference point lies along an existing or 
intended borehole . 

The integrating step preferably further comprises setting 
the predetermined radius around the reference point to the 
radius of the borehole along which the reference point lies. 

This provides the advantage that the integral obtained 
then defines the probability of intersection of two boreholes. 

The method is preferably applied iteratively to a set of 
boreholes and the integrals obtained from the integrating step 
are summed . 

In this way, the summation defines the probability of any 
of the boreholes being within a specified radius of the 
reference point. 

The method may further comprise the step of geometrically 
transforming the second variance matrices, prior to the 
determination of the normalised distribution, into an arbitrary 
co-ordinate system and said integrating step comprises 
restricting the integral to a single dimension so that it 
represents the probability that the borehole lies within a 
specific distance and direction of the reference point. 

In a preferred embodiment, the reference point is the 
position of the bit of a borehole drilling apparatus located 
within a borehole, and the arbitrary co-ordinate system is 
defined by the direction and inclination of the borehole at 
that position. 

This provides the advantage that the integral determined 
in the integrating step defines probabilities in either of the 
directions ahead / behind, left / right or above / below the 
bit. 

According to another aspect of the present invention, 
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there is provided a method of borehole drilling comprising a 
method as defined above, and controlling drilling of the 
borehole according to the probability obtained from the 
integrating step. 

* According to a further aspect of the present invention, 

there is provided a method of determining position uncertainty 
in a directional borehole, the method comprising: 

measuring the azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement of successive points from a 
surface location using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made; 

modelling the relationship of the magnitude of both 
systematic and random errors with azimuth and inclination at 
the measuring locations, and its displacement from the origin 
and elapsed time as a weighting function; 

determining a variance matrix representing position 
uncertainty of each measurement location; and 

obtaining from the variance matrix an ellipsoid bounding 
the volume of space containing the probable position at a 
predetermined level of confidence of each measurement position. 

In a preferred embodiment, the step of determining the 
variance matrix accounts for the contribution to the variance 
matrix of the uncertainty in the position of the surface 
location . 

Preferably, the step of determining the variance matrix 
accounts for the contribution to the variance matrix of the 
correlation between individual systematic error terms in 
successive sections of a borehole survey run in multiple 
sections . 

According to a further aspect of the present invention, 
there is provided a method of determining position uncertainty 

* 

in a directional borehole, the method comprising: 

measuring the azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement of successive points from a 
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surface location using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made; 

modelling the relationship of the magnitude of both 
systematic and random errors with azimuth and inclination at 
the measuring locations, and its displacement from the origin 
and elapsed time as a waiting function; 

determining a variance matrix representing position 
uncertainty of each measurement location; and 

obtaining from the variance matrix an ellipse defining an 
area containing the probable position at a predetermined level 
of confidence of each measurement position, or obtaining the 
probability in a predetermined direction. 

The method may further comprise the step of determining 
the minimum number of survey sections and measurement locations 
required to achieve a predetermined level of confidence by 
iterative application of the method defined above. 

A preferred embodiment of the invention will now be 
described, by way of example only and not in any limitative 
sense, with reference to the accompanying drawings in which :- 

Figure 1 is a schematic illustration of two boreholes and 
ellipsoids of position uncertainty surrounding points on the 
bor eho 1 e s ; and 

Figure 2 is a schematic illustration of the probability 
of collision between a pair of directional boreholes. 

Nomenc 1 a tur e 

V The variance matrix at a particular 

station 

COV The covariance matrix at a particular 

station. 

H The matrix used to represent the 

covariance matrix. 
N, E, V North, East and Vertical respectively 
var (,) The variance of positional 

uncertainties . 
x, y, z x, y, z coordinates 

cov (,) The covariance positional 
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D 
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uncertainties . 

The position vector resulting from the 
measuring errors. 

The magnitude of the individual j th 
systematic measuring error. 

0 

The maximum magnitude of the individual 
j th systematic measuring error. 
The magnitude of the j th random 
measuring error at the i th station. 
The weighting function that takes into 
account that individual measuring 
errors may depend upon the actual 
borehole vector or duration of the 
survey . 

The unit vector indicating the 
direction in which the measuring error 
is effective. 

The measured along hole depth at the i th 
station. 

The magnitude of the individual j th 
measuring error for the i th instrument 
used in a survey comprising consecutive 
runs using different survey 
instruments . 

The total along hole depth of the 
survey obtained using a specific 
instrument in a survey comprising 
consecutive runs using different survey 
instruments . 

The total along hole depth of the 
borehole . 

The expected position of a particular 
station. 

The true Position of a particular 
station. 

The error vector. 

The minimum distance between two 
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stations . 
R The borehole radius. 

cr x 2 The variance of the error distribution 

of the x direction. 
cr w The covariance between the error 

xy 

distributions in the x direction and 
the y direction. 

The correlation coefficient given by 



xy 



/cr x 2 /a y 2 



AtrUj The maximum value of the j th individual 

measuring error. 
X The parameter determining the size of 

the ellipsoid. 
X' The parameter used to scale the 

individual measuring errors. 
P The probability. 

The meanings of a 5/ a.j (i) , aj (1) , and aj (2) are described in the 
text . 

When a borehole is surveyed, imprecision in the 
observations leads to the development of a degree of 
uncertainty in the computed position at each point (or station) 
of the survey. This uncertainty is typically expressed as a 
variance matrix, the elements of which define the variances and 
covariances of the error in position in the North, East and 
Vertical directions. To facilitate visualisation of the 
position uncertainty, ellipsoids of the form 

ArV^Ar = X 2 

where V is the variance matrix and Ar is the position vector 
resulting from the error vectors, which describe volumes of 
equal probability density, may be extracted from the variance 
matrix. 

The basic objective of the improved method is to 
determine the positional uncertainty of each station in a 
borehole survey. This may be represented by either an 
ellipsoidal volume of uncertainty centred about the expected 
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position of the station or by a covariance matrix. The 
covariance matrix at a particular station is determined by- 
taking the expected value of the matrix formed by 
premultiplying the column vector of the positional error about 
the expected value of the error at that station with its row 
vector (i.e. its transpose). In other words: 
COV = E { [Ar - E (Ar)] [Ar - E (Ar)] T } 

This matrix contains a statistical description of the 
probability distribution of the position vector of the station 
under consideration and has dimensions given by: 



COV = 



var (AE, AE) cov(AE,AN) cov(AE,AV) 
cov(AE,AN) var (AN, AN) cov(AN,AV) 
cov(AE.AV) cov(AN,AV) var(AV f AV) 



The eigenvalues and eigenvectors of the covariance matrix 
determine the axes of an ellipsoid containing the possible 
positions of the station when all the individual measuring 
errors are taken into account . 

In the analysis below it is assumed that the probability 
distribution of each source of measurement error is bounded by 
finite upper and lower limits. This assumption implies that 
the ellipsoid will contain all the positions that are possible 
for the station. 

If, on the other hand, it is assumed that the 
distribution of each error is not bounded then, by the Central 
Limited Theorem, the ellipsoid describes the volume containing 
a three-variate normal distribution (i.e. normal distributions 
in each of East, North and Vertical directions) of possible 
station locations . 

Present art leans towards extracting the ellipsoid from 
the covariance matrix at an early stage, thereafter 
manipulating it as a geometrical volume rather than a 
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statistical entity. 



This approach has developed because of 



the difficulty of visualising the abstractions embodied in the 
covariance matrix and has led to a number of difficulties, it 
is meaningless, for example, to .combine error ellipsoids as 
summation of geometric characteristics such as the length of 
axes or the enclosed volume does not reflect the correct manner 
of combining variances of covariances . It is, therefore, 
desirable to retain the covariance matrix* for as long as 
possible; this approach is taken in the methods described 
below. 

Systematic and random measurements errors may be combined 
to produce a covariance matrix which fully describes the 
positional uncertainty of a station. 

If an East, North, Vertical co-ordinate system is assumed 
the positional uncertainty Ar at the first station is given by: 



where the first term on the RHS represents the contribution to 
the positional uncertainty of the systematic error terms and 
the second term represents the contribution of the random error 
terms . 

The individual systematic error terms, Anij, are 
considered to be the same at each station. Their actual values 
for any specific survey are unknown, but are drawn randomly 
from distributions centred on their expected values. The 
values for the individual random error terms, Am^ (1) , vary from 
station to station; these too are drawn randomly from 
distributions centred on their expected values. 

Since, the systematic error terms at the second station 
are the same as the first, the positional uncertainty of the 
second station will be: 



D 



D 
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Ar = ^ Am. j w. (s) v j (s) ds + ]^ A ™J J w j < s ) v j < s > * 

i o i o 

Hence, it follows that the positional . uncertainty at the nth 
station is: 




As mentioned above the covariance matrix at a station with a 
position vector, r, has the form: 

COV - E { [Ar - E (Ar) ] [Ar - E ( Ar ) ] T } 

This may be expanded to give: 

COV = E { ArAr T } - E {Ar} E { Ar T } 

This represents a matrix of the form E {ArAr T } whose centre is 
about the position vector of the station, r, plus the expected 
value of the positional error E { Ar } at r. 
Thus, writing: 



Da 



= Jw j (s)v.(s)ds 



Di 



and 



j ■ 



= J w j (s)v j (s)ds 



Di-1 
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the covariance matrix: 
COV = E{ArAr T } 

on substituting (1) becomes: 

COV =E {^U } Am.a. + X.Aml l \l 1) + ... + ^ j A^ n \l ay ) 

(^Am^^Am/V" + ... + X j Ai/ n) a. (n, ) r } 



which gives : 



COV s= 



>V = Z.Z k E (Am.Ang a.a* + ZjZ^ Am.Am^ Ja.a (I)T 
'j£ k E[ Amj } Amj.Ja^ a£ + . . . 



+ . 



lAm. Am^Jaj 



a k + ... + Z.2 k E 



(aiti.* 



(n) - (n) j (n) 

Am. ya. 



Cn)T 



Since the error terms Am^ {1) . Am j <2) . .., Amj (n> are random, 
these terms must be uncorrelated with each other. If it is 
also assumed that the individual systematic error terms are 
uncorrelated both with each other and with the random terms, it 
follows that the covariance matrix will be: 



COV = L J a j a;E(Am;) + Z ja / 15 a j (l,T E(Am; ,)2 ) 



(2) 



_ CO 
+ ... +£jaj a a 
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The above assumption is valid provided no individual 
measurement of a tool is dependent on any other individual 
measurement . 

Equation (2) represents the covariance matrix formed by 
summing the systematic errors together over the entire length 
of the borehole plus the sum of the covariances matrices formed 
by summing the random errors at each station. The position 
error at any station is, therefore, dependent not only upon the 
position vector of the station, (cf . Wolff and de Wardt) but 
also upon the number of stations before it in the traverse (cf . 
Walstrom et al . ) . The combination of these hitherto 
incompatible approaches is at significant improvement over the 
prior art . 

As the number of stations in a survey increases, so there 
is a greater tendency for random measurement errors to cancel 
one another; the total contribution of random measurement 
errors to the covariance matrix therefore tends to a limiting 
value. This has implications for the planning of borehole 
surveys in that both the optimum station interval and, as will 
be seen from the next section, the optimum number of survey 
sections must be considered. 

There are also implications for borehole surveying 
practices at shallow depths which remained unexposed by present 
art. Current, Wolff and de Wardt based, models not only 
underestimate uncertainty because the contribution of random 
errors is not considered, but also, as the number of stations 
is low, are in error by a changing degree. 

It is a common practice to survey boreholes in multiple 
sections which are then concatenated ("tied" ) into one survey. 
The individual surveys may be taken with tools of the same type 
(but perhaps not exactly the same instrument) , but frequently 
different types of instrument are used in the different 
sections. It is now shown how the results of two or more such 
survey sections may be tied together. The technique may also 
be used to account for uncertainty in the surface location of 
the borehole . 
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Consider first the case where the survey comprises the 
results of only two tools, tied together at a "tie point" of 
depth D. In this case, the position error vector at a station 
at depth D M is given by: 

Ar m XjAinJ? J° w j (s) ^ (s) ds + ^Anif 7 * J ** w. (s) v. (s) ds 

+ random terms 

where the first term on the RHS of the equation represents the 
contribution from the first tool and the second term represents 
the contribution from the second tool . The superscripts (1) and 
{2) refer to the first and second tools respectively. 
If this is rewritten as: 

Ar = E j Am j (1) a j (1) + EjAm-j (2, aj <2) + the random terms 



the covariance matrix: 



COV = E { ArAr T } 



becomes : 



COV = E {Z^a/V^Am^Am 



(i) 

k 



+ ZfaW up Amf Am™ + ^.W z^^Am? 



+ the contribution of the random errors to COV 

since the individual random errors are uncorrelated both with 
each other and the systematic errors. 

It can be seen that the contribution to COV of the 
systematic errors is different depending on the assumptions 
made concerning the error data. If, as in the previous 
section, it is supposed that the individual sources of 
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measuring error for each tool are uncorrelated and ignoring, 
for the purposes of clarity only, the contribution of the 
random errors, the covariance matrix will become: 



COV = Z.a. (I) a<" T 

J J * 



E(Am (,)2 ) + Eja/Haf^Amj « Am, < 2 >) 
+ X.a/ 2 > a^XAm. (2) Am, + Z.af > Am. < 2 >) 



(3) 



If, in addition, it is supposed that each source of measuring 
error in the first tool is uncorrelated with the corresponding 
error of the second tool, then the non- leading terms become 
zero and COV becomes : 

COV = Z.a/ , V (,jT E(Am (,)2 ) + 2.a/«af T E (Am/ 2 > 2 ) 

Extending equation (3) to consider a survey comprising n tied 
traverses, the equivalent form is: 

♦vr'j" 1 ^ 4 ","'^/")*-.. 

+ ... 



(4) 



Note that the ratio between the number of non-leading 



where 



terms to the number of leading terms, E^a^ <p) a.j <p)T E 
1 < p < n, increases as n increases. This implies that the 
discrepancy between the correlated and uncorrelated forms of 
the covariance matrix will increase as n increases and hence it 
is important to establish the relationship between the error 
terms of the different tools when the covariance matrix for a 
tied survey is computed. Note also that, since the magnitude 
of the vector (p> is dependent upon the length of the section, 
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p, so too is its contribution to the covariance matrix. 

Also seen in (4) is an implication for the design of 
borehole surveys which is not exposed by present art, namely 
that by careful selection of the type of survey instrumentation 
such that the uncorrelated form of the matrix is appropriate, 
then dividing the borehole into a number of length sections , 
each making approximately the same contribution to the 
covariance matrix, will decrease the total error at depth. 

Ellipsoids (or in two dimensions ellipses) of equal 
probability density can be denoted in vector notation by: 

Ar T COVAr = X 2 

According to the method of Wolff and de Wardt hitherto 
employed a similar relationship exists between the variance, 
a 2 j7 and the quoted possible maximum error magnitudes, AMj , of 
each measurement error 

(AMj) 2 = X' 2 ^ 2 

They proceed, erroneously, to equate X and X' to produce a new 
matrix, H, such that 

A^H^Ar = 1 

The eigenvalues of H are the components u, v and w in each of 
three, orthogonal directions of the maximum position error. 
The maximum position error itself is clearly 

Ar^ = /u 2 + v 2 + w 2 

which does not lie within the ellipsoid the axes of which are 
u, v and w. 

The approach of Wolff and de Wardt is therefore rejected 
in favour of a new method, described below, which places the 
extraction of ellipsoids or ellipses on a firmer statistical 
basis and, by establishing the values of X and X' directly, may 
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be used to extract ellipsoids and ellipses of any desired 
confidence level . 

Unlike the method of Wolff and de Wardt, in which the 
measuring errors are quoted as maximum possible errors, it is 
assumed that they are normally distributed around a mean of 
zero. This assumption is expected and borne out by 
experimental evidence. Even if not generally provable, the 
combination of error distribution will, by the Central Limited 
Theorem, tend to produce a' normal distribution for the 
resulting position error. Furthermore, the assumption 
acknowledges the fact that it is difficult to ascertain exactly 
wha; the maximum individual measuring errors are. 

The two definitions, although differing, may be related. 
Suppose the jth measuring error is thought to lie between 
limits of ±AM j . Since Am-j is assumed to follow a normal 
distribution with mean zero it implies that k = AM^/a (which is 
equivalent to X') is a standard normal statistic. Thus, if it 
is desired to establish that Amj lies between the desired limit 
upto a high degree of confidence it is simply a matter of 
choosing a value of k from the standard normal tables which 
gives a high probability. For example k = 4 gives a 
probability of close to 99.99%. 

Now, if each individual measuring error is normal, it 
follows that the region of uncertainty associated with the 
measurement of a point will follow a 3-variate normal 
distribution. Suppose, more generally, that 

4 

X = ( X-^ t X 2 , . . . Xp ) 



is a co-ordinate in p-space, which follows a p-variate normal 
distribution with mean \i and covariance matrix COV, i.e. x ~ JSL 
(/i, COV) , and where each individual component Xj (1 <. i< p) 
follows a normal distribution, i.e. Xj ~ N 1 (ij,cr 2 i • 

The multi-variate normal distribution will have constant 
probability densities on surface of the form:- 



(x - /x) T COV" 1 (x - /z) = X 2 



(5) 
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It can be shown that the U statistic 
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U = (x - /x) T COV" 1 (x - fi) 

follows a Chi squared distribution of p degrees of freedom, 
i.e. U - x P 2 - Thus the form of a surface describing a region 
of uncertainty in p dimensions can be determined upto any- 
specified confidence level using (5) where X; is given by 

which in three dimensions simplifies to 

X = / X 2 3 

Clearly, the dimensions of the ellipsoid of desired probability 
are therefore the square root of the eigenvalues of X 2 COV. 

The ellipse defining the boundary at a given confidence 
level in the horizonal (XY) and vertical (XZ and YZ) planes may 
be obtained by extracting the appropriate 2 by 2 matrices from 
COV and using the Chi square statistic for 2 degrees of 
freedom. Likewise the confidence level along any of the X, Y 
and Z axes may be determined using the Chi square statistic for 
one degree of freedom. 

To determined the ellipse of uncertainty in an arbitrary 
plane, i.e. one other than XY, XZ, YZ, it is necessary to 
integrate the probability density function described by the 
covariance matrix with respect to the direction normal to the 
plane of interest. The ellipse thus produced is, it should be 
noted, neither a section nor a projection of the ellipsoid. 
Although such a method can be implemented, the same result can 
be obtained by transforming COV into a new co-ordinate system, 
X'Y'Z', in which one of the axes in aligned in the desired 
direction. The eigenvalues and eigenvectors of the matrix thus 
obtained may then be used in the manner described above. The 
validity of this approach can be seen by considering the case 
where the desired co-ordinate system is that of the ellipsoid 
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itself. In this circumstance, COV is transformed into a 
diagonal matrix, the leading terms of which are the variances 
in the directions corresponding to the axes of the ellipsoid 
and the eigenvectors are X' , Y' and Z' . Since such 
transformations are orthogonal the eigenvalues remain, of 
course, those of COV. 

Referring to Figure 1, points r x and r 2 are points along 
boreholes A and B respectively, which have been acquired using 
some borehole surveying technique. 

The true position of well A, r[, will lie somewhere in 
the ellipsoid centred about r x and the true position of well B, 
r^, will lie somewhere in the ellipsoid centred about r 2 . Note 
that the centres of the ellipsoids will be offset from the 
boreholes if any biased systematic terms, such as a drill 
string magnetization error, are considered in the determination 
of the position uncertainty. If the radii of the boreholes are 
Ri and R 2 respectively then, since the direction of the well 
paths are not known precisely, all that can be said about the 
surfaces of the boreholes perpendicular to the direction of the 
well paths at these points is that they lie somewhere in 
spheres of radius R x centred about r x and R 2 centred about r 2 
respectively. Hence it is certain that the boreholes will not 
coincide at these particular points provided d = |T| , the 
distance between r x and r 2/ is greater than the sum of the 
borehole radii . 

Assuming that the individual sources of measuring error 
for a survey tool are approximately normally distributed and 
independent of one another, which assumption is borne out by 
practical observation, it follows by the Central Limit 
Theorem, that each of the coordinates at r^ = (x L/ y lf z[) will 
tend towards being normally distributed with expected values 
~i/ 7n ~i* and variances a xl 2 , a yl 2 , a zl 2 , (say) . Similarly each 
of the coordinates at r 2 = (x 2 , y^, "z 2 ) will tend towards being 
normally distributed. 

Thus, if it is assumed that the points r x and ~ 2 are 
independent of one another then by the properties of combining 
normal distributions it follows that the co-ordinates of the 



WO 96/35859 



PCT/GB96/01123 



22 

vector Y = (x, y, z") given by r 2 - r ± (i.e. the position vector 
which joins the two points at and r* 2 ) will each be normally 
distributed with expected values of x, y and z respectively and 
variances a x 2 = a xl 2 + a^ 2 , o 2 = a yl 2 + c y 2 and o 2 = o z 2 + a z2 2 

Now the expected value E(r) of the random variable 

r = r^ - rj [ = (x, y, z) T 

is 

E(r) = (E(x) ,E(y) ,E(z) ) T = (x 2 - x 1# y 2 V Yl , z 2 - Zl ) T 

and the variance matrix (or covariance matrix) , V(r) , is 

E(r) = E{ [r-E(r)] [r-E(r) ] T } 
= E { - E (r 1 )] [r 1 -E(r 1 )] T } + E { [r 2 -E (r 2 ) ] [r 2 -E (T 2 ) ] T } 





+ 


2 

^x2 ^xiyi + ^X2y2 ^xi zi 


^x2z2 


^xlyl 


+ 


^x2y2 ^yl + ^y2 ^ylzl + 


V y2z2 


^xlzl 


+ 


2 

°x2z2 Gylzl + °y2z2 °zl + 





if r^ and r 2 are independent 

Thus, given that r = (x, y, z) T where x, y and z follow 
the normal distributions 



x ~ N(x 2 - x x , a xl 2 + cr x2 2 ) 



y ~ N(y 2 - y lt o y 2 + cr y2 2 ) 



z ~ N(z 2 - z lt a zl 2 + <x z2 2 ) 

it is assumed that "r follows a 3-variate normal 
distribution with mean E(r) and variance V(r) 

The distribution of d = |r| = Vx 2 + "y 2 + z 2 (i.e. the 
distance between the two points) can be evaluated by taking the 
3-variate normal distribution of r and making the substitution 
x = d cosdslncp , y - dcosfl cos <f>, z = dsinfl . Thus, it follows 
that the probability that the boreholes will not coincide at 
these points is at least as much as 
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P (borehole will not coincide) = 



|° f J f(d,e,(I))d 2 cos8 dddecty 



where 



f (x f y r z) = 1 . exp {- a . r'VtrJ^r} 

(2?r) 3/2 . (detV(r))* 

Referring to Figure 2, it can therefore be seen that the 
probability of collision can be estimated by 

P (collision) = 1-P (boreholes will not coincide) . 

The position uncertainty at any given point along the 
reference well ahead of the current depth can considered to be 
due to contributions due to errors in the survey down to the 
current depth plus those which accumulate over the interval to 
the point of interest. The precision with which azimuth and 
inclination may be set and maintained is analogous to the 
heading reference and inclination systematic error terms of a 
survey tool error specification. Likewise, the circumstance 
considered is similar to that where a wellbore has been 
surveyed in two sections, each using a different survey tool. 
The method described above for computing the position 
uncertainty of such a survey is directly applicable. Note that 
the sources of error in the two sections are completely 
unrelated and so the uncorrelated form combination is the more 
appropriate . 

One method of determining sensible values for drilling 
precision is by analysis of a database of expected and achieved 
results, suitably classified by operation, BHA type, depth, 
formation etc. Another is to consider the variances in azimuth 
and inclination observed when suitable mechanical simulations 
of drilling assemblies are compared against actual results. 
Such a process is typically used in the calibration of 
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mechanical simulations and is often known as " history 
matching 11 . 

In practical circumstances the probability of 
intersection of one borehole with another is "of limited value, 
rather it is the probability of intersection with any of the 
boreholes beneath a drilling facility which is of interest. 
This may be computed according to the following method. 

Let Pi, p 2 # • - • -P n be the probability of a point on the 
reference borehole being within a specified radius of the 
points of closest approach to that point of object wellbores 
l...n< i.e. the probability of intersection with each of the 
object wellbores. The probability of not intersecting the 
first object wellbore is simply p 1 1 = l-p 17 that of not 
intersecting the second object wellbore p^ = 1 - p 2 etc. The 
probability of not intersecting any of them, P no intersectiori , is 
thus 

P- p 1 p 2 p 1 
no intersection "~ r l r l ' " • r n 

and the probability of intersection with at least one, 

^intersection' 

^intersection = l"~^no intersection 

The method described above allows the determination of 
the absolute probability of intersection of a borehole with 
those surrounding it . For the purposes of well planning and 
directional control, however, it is desirable to identify 
positions of least risk of intersection and to navigate through 
those positions, such positions not necessarily being 
equidistant between the existing boreholes. 

Considering any arbitrary point the method described can 
be applied by reducing its position uncertainty to nothing and 
reducing the tolerance around it to zero. This may be achieved 
assigning an empty variance matrix, i.e. a matrix with all 
terms zero, to the reference point and setting R x to zero. The 
probability of the closest points on each of the set of object 
wellbores being within the tolerance R 2 of this point may then 
be computed and summed as shown in the proceeding section. 

By altering the value of R x to any appropriate number the 
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method described above may be used to evaluate the probability 
that a borehole will lie within a given tolerance of any fixed 
position, for instance a geological target or similar feature. 

Whilst drilling it is often desirable to know the 
likelihood that the current reference borehole position lies 
along a specific direction with respect to the object borehole 
or, indeed, some other point such as the intended target. For 
example, it may be important for the reference borehole to pass 
beneath, or to the East or to the North, of the object well. 

Consider the random variable z = z 2 - z x which is the Z- 
co-ordinate of the actual position vector between point A on 
the reference borehole and point B on the object borehole. If 
Z is taken to represent vertically downwards then z > 0 implies 
that the vertical depth of B is greater than that of A. Hence > 
if it is assumed that the vector [x, y, z] T follows a 3-variate 
normal distribut ion the probability that B lies below A is 
equivalent to the probability that z > 0 which can be 
determined as follows 



/*oo to CG 

z>0) = J - m0 J - J 



P(z>0) = J - sQ J - =00 J x=-» f<*/ Y t z) dxdydz 

Similarly it is possible to determine the probability 
that A is North of or East of B by considering the probability 
that y > 0 and z > 0 respectively. 

In fact, one is not limited to this. It is also possible 
to determine what the likelihood of the relative position of 
the object well is in terms of any combination of North, East 
and Above . For example , it can be determined that the 
reference borehole is North and East of the reference borehole 
or even North, East and Above. 

Neither is one limited in terms of North, East and Above. 
For example, when reference and object wellbores approach one 
another normally it is useful to determine the relative risks 
of passing to the left or right. Such a circumstance may occur 
when drilling an horizontal well in proximity to a vertical 
exploration or delineation well. By defining a new frame of 
reference (X' , Y' , Z' ) such that Z' is aligned with the 
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reference borehole vector, X' is normal to this vector in the 
horizontal plane and Y' is normal in the vertical plane then a 
matrix, M, can be derived which defines the transformation 
between (X', Y' , Z' ) and (N,E,V). Applying this transformation 
to the variance matrix such that 
V 1 = M T VM 

produces a transformed variance matrix which, when the method 
defined above is applied, will produce probabilities of "ahead 
/ behind 1 ', "left / right 11 or "up / down" or any combination 
thereof . 

A property of V 1 is that its eigenvalues will be the same 
as those of V although its eigen vectors will, of course, 
differ. Further, if the orientation of (X', Y' , Z') is chosen 
such that they align with the eigenvectors of V then V 1 will be 
a diagonal matrix, the leading terms of which are the variances 
in the directions corresponding to the axes of the ellipsoid of 
uncertainty and the eigenvectors X' , Y' and Z' . 

Using the method described above, charts, maps, strip 
charts or engineering drawings showing the position of a 
borehole or proposed borehole, known as the "reference well", 
and a suitably selected group of other boreholes, such 
boreholes being known as "object wells 1 ', are created. Criteria 
for membership of the set of object wells may include, amongst 
others, being drilled from a specific structure, lying within 
a specific field or coming within a specific distance of the 
reference well . The drawing may display the well in either 
plan or vertical or isometric or other projection. The drawing 
may also display, instead of the actual trajectory of the 
boreholes, their relative separation; such charges have been 
called "travelling cylinder" or "rose diagram" plots. The 
drawing may display ellipses or ellipsoids of uncertainty at a 
specified level of confidence. The probability of intersection 
with each object well may be displayed as an attribute of the 
object well or the joint probability intersection with any of 
the object wells as an attribute of the reference well. The 
choice of display and the attribute employed for representation 
may best be decided according to the task at hand. Examples 
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include varying colour, weight or style of the line 
representing a well; using "wiggle traces" (either plain or 
filled) or more complex denotations such as spot sizes / colour 
or "tadpoles". Typically such drawings are rendered upon 
paper or other tangible and reproducible medium although 
transient displays upon a visual display unit may also be used. 

The evaluation of the variance matrices and the 
computation therefrom of the probability of intersection and 
the analysis and interpretation thereof are generally performed 
numerically by means of a computer program. To be of practical 
use such a program should also encompass a. means of storage and 
retrieval of a large number of borehole surveys, together with 
a means of computing the co-ordinates of stations therein (at 
least according to the minimum curvature method) and of 
interpolating the position and borehole vector of points 
between stations. It should also provide a facility to create 
and maintain associations between borehole surveys and error 
models. The error models should contain at least a description 
of the source and magnitude of the measurement errors, nij and 
m j U) , applicable to the type of instrumentation used to gather 
the data. In the circumstance where the survey is acquired in 
two or more sections such associations mst be maintained for 
each section of the survey. 

In the circumstance that large numbers of survey are 
stored, examples being large offshore platforms or Arctic 
drilling locations, the program should provide a facility for 
the classification of boreholes and their surveys and their 
extraction from storage according to the classification. Such 
a facility simplifies the creation of sets of object wells and 
minimises the risk that relevant data will be inadvertently 
neglected. 

Furthermore, the program should provide a facility for 
the creation of, at least, a description of the trajectory for 
planned boreholes in a manner analogous to that used for 
surveys of actual wells. It should also provide for division 
of the trajectory of a planned borehole into survey sections 
and for the association of each section with an error model . 
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In practice it is useful to allow the creation of many such 
divisions and associations in order to simulate, compare and 
contrast the influence of possible surveying programmes upon 
borehole position uncertainty at different stages of the 
drilling process. 

To support this function it is desirable that a means is 
provided to specify the surface location of drill sites and the 
position and orientation of specific wellheads (slots) therein 
and the uncertainty in position of the drill site and the 
slots. 

In addition it is necessary for such, a program to provide 
a function for the computation of the separation of boreholes 
one from another in the horizontal plane, the plane normal to 
the borehole vector and three dimensions and the determination 
therefrom, by iterative or other means, of the points of 
closest approach of the reference well with each member of the 
set of objects wells. 

Each of these attributes of a practical program may be 
implemented in a variety of ways using a variety of algorithms. 
They may be implemented within a single program or as a suite 
of programs sharing common data. In a preferred implementation 
the features are provided as a suite of inter-related programs 
sharing a common data and derived from common subroutine 
libraries . 

In a preferred implementation the set of errors rtij and 
m/ 11 for magnetic survey tools are identified as: 

1. Error in magnetic declination 

2 . Magnetic compass error 

3 . Inclinometer error 

4. Misalignment with the centre line of the borehole 

5 . Along hole depth measurement error 

The weighting function for error term (2) above accounts 
for the sensitivity of methods of compensation for drill string 
magnetic interference to uncertainty in determination of the 
Earth's magnetic field strength and dip angle. 

For gyroscopic tools the errors are identified as: 
1. Error in alignment with reference direction 
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2 . Gyro compass error 
3 • Inclinometer error 

4 . Misalignment with the centre line of the borehole 

5 . Along hole depth measurement error 

The weighting function for error term (2) above accounts 

for the time dependency on these errors in "continuous" survey 

tools. 

For inertial tools the errors are identified as: 

1. Error in alignment with reference direction 

2. X, Y and Z displacement errors 

The weighting function for error term (2) accounts for 

the time dependency of these errors . 

For all tools the following errors are identified: 

1. Displacement of tool from centre line of the 

borehole . 

The errors enumerated are representative only as the 
method presented can account for an arbitrary set of errors . 

An extended version of this preferred implementation 
incorporates the facility to describe drilling precision in a 
manner suitable for the determination of the wellbore position 
uncertainty. 

A specific method of using a typical least risk of 
borehole intersection chart created according to the method 
presented is given below. Such a chart may be used in the 
initial stages of well planning to identify the safest 
positions to which to navigate. 

1. Create base charts or maps of the location and, 
optionally, the associated position uncertainty of the set of 
object wells at specific vertical depths. 

2. The area covered by the base map, or any portion 
thereof, is divided into cells, the shape and size of which are 
appropriate to the scale of the base map and appropriate to the 
task in hand. Each cell is aligned a value determined as the 
probability of borehole intersection of its centre. 

3 . The surface thus defined may be rendered according 
to any appropriate method including, but not limited to, 
contouring, shading or colouring. 
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4. By examination, either visual or otherwise (for 
example by use of algorithms to determine surface minima) , of 
a set of such charts created at regular depth intervals the 
optimum position of each depth for a new reference well may be 
determined. 

5, A modification of the above method is appropriate 
for visual displays wherein the base map of step (1) and the 
cells defined in step (2) are three dimensional and rendered 
using appropriate 3 dimensional visualisation software. 

A process of drilling incorporating the new method is as 
follows : 

1. Determine the desired course of the borehole 
according to the procedure described above or any other 
appropriate to the circumstance at hand. 

2 . Set forth a safe limit on the allowable risk of 
borehole intersection according to either the relevant health, 
safety or environment regulations extant at the time and place 
of drilling or the financial risk, which ever is the lesser. 

3. Drill an initial length of the borehole, the BHA 
and drilling parameters (particularly weight-on-bit and rotary 
speed) used being selected to give adequate control for the 
section of well. 

4 . At regular intervals along the borehole take 
observations of the borehole inclination and azimuth using 
either magnetic or gyroscopic instruments and hence determine 
both the position of the bit and its associated position 
uncertainty, 

5. Using a mechanical simulation of the BHA or 
otherwise estimate the trajectory of the borehole over the next 
length to be drilled and determine first the position and 
associated position uncertainty of the drill bit at the end of 
this length according to the drilling tolerances of the BHA and 
the quality of the estimation and thence the risk of 
intersection. 

6. If the risk of intersection exceeds the safe limit 
set forth in step (2) , repeat step (5) successively considering 
alternative trajectories achievable with the BHA in use and 
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select the safest trajectory compatible with reaching the 
desired objective and remaining within other drilling 
constraints, for example torque or drag limits. 

7. If the safest achievable trajectory exceeds the 
safe limit, repeat step (5) successively considering 
alternative BHA configurations. 

8. If the safest achievable trajectory exceeds the 
safe limit, repeat step (7) considering trajectories achievable 
from prior points along the borehole. 

A preferred embodiment of this process offers advantage 
over prior art by permitting the computed risk of intersection 
to be associated with cost related data, for example the 
operating cost per hour of the drilling rig. Using this 
preferred embodiment decisions such as whether or not to 
extract and replace the drilling assembly or the depth at which 
to plug and sidetrack a borehole the trajectory of which cannot 
be adequately controlled, both of which are expensive 
operations, can be made on an objective, quantative basis ; 
presently such decisions are made by subjective judgement which 
may, retrospectively, prove incorrect. 

It will be appreciated by persons skilled in the art that 
various alterations and modifications are possible without 
departure from the scope of the invention as defined by - the 
appended claims. 
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CLAIMS : 

1. A method of determining the risk of a collision of 
a borehole, the method comprising: 

measuring azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement from a surface location of 
successive points using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made and their relationship with 
azimuth and inclination of the borehole and displacement from 
the origin and elapsed time at each successive point of 
observation ; 

determining a first variance matrix representing position 
uncertainty at each location at which said measurements are 
made ; ' 

determining the point of closest approach of the borehole 
to a reference point; 

determining second variance matrices representing the 
position uncertainty at said reference point and said point of 
closest approach; 

determining a vector between the reference point and the 
point of closest approach; 

determining a normalised distribution of the length of 
said vector according to the second variance matrices; and 

integrating the normalised distribution between limits of 
zero and a value equal to the sum of a predetermined radius and 
the radius of the borehole, wherein the integral obtained 
represents the probability that the borehole lies within the 
predetermined radius of the reference point. 

2. A method according to claim 1, wherein the step of 
determining the first and / or second variance matrix and / or 
the normalised distribution are performed by application of a 
stochastic technique, and preferably by application of Monte 
Carlo modelling, and the integrating step is performed by 
examination of the frequencies thus derived. 

3. A method according to claim 1 or 2, wherein the 
first variance matrix is determined from estimates of expected 
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drilling precision. 

4 . A method according to any one of the preceding 
claims, wherein the reference point lies along an existing or 
intended borehole. 

5 . A method according to claim 4 , wherein the 
integrating step further comprises setting the predetermined 
radius around the reference point to the. radius of the borehole 
along which the reference point lies. 

6 . A method according to any one of the preceding 
claims, wherein the method is applied iteratively to a set of 
boreholes and the integrals obtained from the integrating step 
are summed . 

7. A method according to claim 6, further comprising 
the step of geometrically transforming the second variance 
matrices, prior to the determination of the normalised 
distribution, into an arbitrary co-ordinate system and wherein 
said integrating step comprises restricting the integral to a 
single dimension so that it represents the probability that the 
borehole lies within a specific distance and direction of the 
reference point . 

8 . A method according to any one of the preceding 
claims, wherein the reference point is the position of the bit 
of a borehole drilling apparatus located within a borehole, and 
the arbitrary co-ordinate system is defined by the direction 
and inclination of the borehole at that position. 

9 . A method of borehole drilling comprising a method 
according to any one of the preceding claims, and controlling 
drilling of the borehole according to the probability obtained 
from the integrating step. 

10. A method of determining position uncertainty in a 
directional borehole, the method comprising: 

measuring the azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement of successive points from a 
surface location using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made; 
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modelling the relationship of the magnitude of both 
systematic and random errors with azimuth and inclination at 
the measuring locations, and its displacement from the origin 
and elapsed time as a weighting function; 

determining a variance matrix representing position 
uncertainty of each measurement location; and 

obtaining from the variance matrix an ellipsoid bounding 
the volume of space containing the probable position at a 
predetermined level of confidence of each measurement position. 

11. A method according to claim 10 , wherein the step of 
determining the variance matrix accounts for the contribution 
to the variance matrix of the uncertainty in the position of 
the surface location. 

12. A method according to claim 10 or 11, wherein the 
step of determining the variance matrix accounts for the 
contribution to the variance matrix of the correlation between 
individual systematic error terms in successive sections of a 
borehole survey run in multiple sections . 

13 . A method of determining position uncertainty in a 
directional borehole, the method comprising: 

measuring the azimuth and inclination of the borehole at 
intervals along it using magnetic or gyroscopic instruments, or 
measuring the total displacement of successive points from a 
surface location using an inertial instrument; 

determining the magnitude of systematic and random errors 
inherent in the measurements made; 

modelling the relationship of the magnitude of both 
systematic and random errors with azimuth and inclination at 
the measuring locations, and its displacement from the origin 
and elapsed time as a waiting function; 

determining a variance matrix representing position 
uncertainty of each measurement location; and 

obtaining from the variance matrix an ellipse defining an 
area containing the probable position at a predetermined level 
of confidence of each measurement position, or obtaining the 
probability in a predetermined direction. 
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14. A method according to any one of claims 10 to 13, 
further comprising the step of determining the minimum number 
of survey sections and measurement locations required to 
achieve a predetermined level of confidence by iterative 
application of the method. 
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